How to measure a spinon Fermi surface 
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We propose an experiment to identify the potential existence of a spinon Fermi surface by looking 
for oscillatory coupling between two ferromagnets via a spin liquid spacer. Three candidate spin 
liquids are investigated, and it is found that in all cases, long period oscillations should be present, 
the period of which would identify the Fermi wavevector of the spinon surface. 
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In 1973, Anderson [1] proposed the possibility of a spin 
liquid, a state where long range magnetic order is sup- 
pressed by frustration, low dimensionality and/or quan- 
tum fluctuations. In 1987, he resurrected this idea in 
the context of high temperature cuprate superconduc- 
tors [2]. The proposed ground state was a so-called uni- 
form resonating valence bond (RVB) state which pos- 
sesses a Fermi surface for spin excitations. Although it 
turns out that the undoped cuprates are antifcrromag- 
nets, the spin liquid concept in the context of doped 
cuprates is a very active field of study [3] . 

Frustration plays a major role in suppressing mag- 
netic order, and the original RVB idea was devel- 
oped for a triangular lattice. Since then, other frus- 
trated lattices have been identified, including pyrochlore, 
kagome, and hyper- kagome [4]. In the past few years, 
three candidate 5 = 1/2 spin liquids have been identi- 
fied [5]: K-(BEDT-TTF)2Cu2(CN)3 (distorted triangular 
lattice) [6], ZnCu3(OH)6Cl2 (kagome) [7], and Na4lr308 
(hyper-kagome) [8]. All of these materials are insulat- 
ing, have a large Curie- Weiss temperature, and yet show 
no ordering down to the lowest temperatures measured. 
In each case [9-13], a uniform RVB state has been pro- 
posed as the ground state with a resulting spinon Fermi 
surface. There is indirect evidence for such a surface. 
The BEDT salt exhibits a linear T specific heat [14] 
characteristic of a Fermi surface. The herbertsmithite 
ZnCu3(OII)6Cl2 exhibits a sub-linear T dependence of 
the specific heat [15], but this can be understood from 
self-energy corrections [9]. In addition, a Pauli-like sus- 
ceptibility is inferred at low temperatures once the effect 
of impurity spins has been factored out. 

On the other hand, this evidence for a spinon surface 
is indirect. Spin glass behavior, or a magnon dispersion, 
w^, with a power z equal to the dimensionality, can also 
lead to a linear T specific heat [4]. Therefore, it would be 
desirable to have a more direct test. Now, the existence 
of a Fermi surface implies the presence of Friedel oscilla- 
tions. For a spin liquid, this would not show up in the 
charge channel, but would show up in the spin channel. 
The challenge is how to detect this. 

As is well known, oscillatory coupling has been seen 
between two ferromagnets separated by a paramagnetic 
spacer (Fig. 1) [16]. This is a consequence of Kohn 
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FIG. 1: (Color online) The proposed experiment involves two 
ferromagnetic layers (F) with a spin liquid spacer of variable 
thickness, z. Depending on the sign of the oscillatory cou- 
pling, the two ferromagnets will be aligned or anti-aligned. 

anomalies of the Fermi surface which appear in the spin 
susceptibility [17], but are difficult to observe by neutron 
scattering because of their weak intensity. These anoma- 
lies are due to extremal spanning vectors of the Fermi sur- 
face. There are several types: 2kp vectors, umklapp vec- 
tors, and vectors which connect different Fermi surfaces. 
Large vectors are difficult to observe in the oscillatory 
coupling since their period is comparable to the lattice 
constant [9] (and are therefore damped by roughness of 
the layers), but small vectors have been prominently ob- 
served in transition metal multi- layers [16] . The existence 
of small vectors (long period oscillations) is generic, as 
large Fermi surfaces have small umklapp vectors, whereas 
small Fermi surfaces have small 2kp vectors. Oscilla- 
tory experiments have the additional advantage of being 
able to detect multiple q vectors with different orienta- 
tions depending on the growth direction of the layers. 
These oscillations are also strongly suppressed by a spin 
gap [18], as the latter leads to an exponential decay of 
the susceptibility with a distance scale set by the lattice 
constant [19]. 

We will now apply this idea to spin liquids. We start 
with the BEDT salt, originally discovered by Geiser et 
al. [20]. As discussed by Shimizu et al. [6], this ma- 
terial is composed of dimers that sit on a distorted 
triangular lattice. This particular salt, though, has a 
very small anisotropy of the hopping integrals of 6%. 
The resulting Fermi surface is shown in Fig. 2a, which 
was derived from the dispersion — 2t' cos(fcb6) -I- 
4tcos(fch6/2)cos(fccc/2) - ^ with i=54.5 meV, i'=57.5 
meV, and ijl—-A6.2 meV to achieve half-filling [b and c 
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FIG. 2: (Color online) Spinon Fermi surface for (a) 
K-(BEDT-TTF)2Cu2(CN)3 (6=8. 59A, c=13.40A) and (b) 
ZnCu3(OH)6Cl2 (a=6.84A). Spanning vectors are indicated 
by arrows. 

are the orthorhombic lattice constants parallel to the lay- 
ers). This surface is removed by a Mott transition [21], 
but the transition is suppressed by pressure. We as- 
sume that the spinon surface in the Mott phase is the 
same as the 'band' Fermi surface [22]. We note that the 
shortest spanning vector along kc is an umklapp vector 
of length 0.947r/c=0.22A~^, which would give rise to a 
real space period of 28. 5A, comparable to that seen in 
multi-layers [16]. We should remark that the real Fermi 
surface (under pressure) has been seen both by angle de- 
pendent magnetoresistance and Shubnikov-de Haas oscil- 
lations [23]. 

We can perform the same exercise for the herbert- 
smithite [7] . This material has copper atoms on a kagome 
lattice, with three copper atoms per unit cell. Assum- 
ing a near-neighbor spinon hopping, t^, one has three 
bands, the middle of which is half-filled. Its disper- 
sion is given by E/tg ~ 1 ~ X — fi, where AX^ — 
*12 +*13 +*23 + 3ti2ii3i23 with tu = 2 cos(fcy a/2) , i 13 = 
2cos(fcya/4 -I- V3fca;a/4),<23 = 2cos(fcya/4 - V^k^a/A) 
and a the lattice constant (each bond has a length a/2). 
The spinon Fermi surface, previously derived by Ma and 
Marston [11], is shown in Fig. 2b. It has a short umklapp 
spanning vector of 0.50(27r/\/3a)=0.265A~^, giving rise 
to a period of 23. 7A, similar to that estimated for the 
BEDT salt. 

One issue with these two materials is that they are 
two dimensional. To see an oscillatory period as dis- 
cussed above requires that the spanning vector have a 
component perpendicular to the multi-layer (i.e., along z 
in Fig. 1). This would require growing the material with a 
surface not perpendicular to the '(001)' direction, which 
may be a difficult undertaking. Along with this would 
be the probable difficulty of growing these materials with 
controlled thickness between the two ferromagnets. 

We now turn to the hyper-kagome case. Na4lr308 has 
several advantages. It is cubic, and is also a transition 
metal oxide. As such, it should not have some of the 
growth difficulties mentioned above. This material has 



12 Ir atoms in the unit cell [8]. Each Ir atom has 4 Ir 
neighbors, forming a network of corner sharing triangles. 
Assuming again near-neighbor hopping, the resulting 12 
bands can easily be found by numerical diagonalization. 
The spinon Fermi surface for this case has been shown 
by Zhou et al. [12], and consists of two hole pockets 
around the R point of the simple cubic zone, and one 
electron pocket around the F point (Fig. 3a). The sur- 
faces are small, and thus the shortest spanning vectors 
are of the 2kp variety. To a first approximation, we find 
that the surfaces can be approximated by spheres (the 
cubic anisotropy is of order 10%). The kp vector for 
the electron surface along (100) is 0.3l7r/a and that of 
the hole pockets 0.267r/a. As there are two hole pockets, 
their contribution will dominate over the electron pocket 
in the oscillatory coupling. With a of 8.985A, this results 
in a predicted oscillatory period of 34. SA. 

We can make a quantitative estimate of the oscillatory 
coupling following the literature in the GMR (giant mag- 
netoresistance) field [17]. We assume that the coupling 
of the spins in the ferromagnet to those in the spin liq- 
uid can be described by an effective contact interaction 
yl(5(r — ri)s- Si where s is the spin vector of the spin liquid 
at position r and Si the spin vector of the ferromagnet 
at position r^. Values for A can be estimated by calcu- 
lating the transmission properties of the barrier between 
the ferromagnet and spacer layers, although these calcu- 
lations are quite involved [18, 24]. This effective contact 
interaction then leads to a coupling between a site in the 
first ferromagnet layer to one in the second 

J., = J rf^qx(q)e- (--) (1) 

where V is the unit cell volume of the spin liquid, and 
x(q) the static susceptibility of the spin liquid in units of 
2/i^ per unit cell. The interlayer coupling per unit area 
is then obtained by summing over the lattice sites in the 
second ferromagnetic layer, 

I^ = P'EJ^^ (2) 

j 

where d is the layer spacing of the spin liquid. Substitut- 
ing for Jij , we have 

""11 

(3) 

with z oriented as in Fig. 1 and the sum is over r vectors 
in the second ferromagnetic layer. By Bloch's theorem, 
the sum reduces q|| to zero, and we obtain 

= I dq.e^^'^xiOAq.) (4) 

To proceed further, we have to specify x- As in the spin 
liquid literature [9] , we will assume that x is given by the 
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bare polarization bubble times the Gutzwiller projection 
factor gj = 4 [25]. This then yields 



/(z) = - 



327r4T/2 



;/(£k) - /(ek+q) 



fk 



(5) 

where q = {0,0, q^), f is the Fermi-Dirac function, ek 
the spinon dispersion, and g the effective g-factor of the 
spin liquid (5 is the magnitude of the spin vector of the 
ferromagnet, and s the magnitude of the spin vector of 
the spin liquid). Note that there is an implicit double 
sum over band indices. 

The integrals over kz and qz can be converted to inte- 
grals over Ek and ek+q, and the remaining integral over k|| 
will be dominated by extremal q vectors connecting the 
spinon Fermi surface(s). The variation of the extremal q 
vector is given by an expansion of the spinon dispersion 
about kp [26] 



ek = Vz{kz 



+ ^{klD^X + k^Dyy + 2k^kyD^y) (6) 



where Vz is the Fermi velocity along z and D is the inverse 
mass tensor. In our cubic case, D^^ — Dyy. Setting Ck 
to zero and solving, we find 



kz — ki 



1 

2^ 



{(kl + kl)D^^ + 2h,kyD^y) (7) 



This is easily diagonalized by rotating 45 deg rees in kx , ky 
space 



kp 



1 
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k^D2 



(8) 



i any- 



where Di2 

Because of the multiple Fermi surfaces in Fig. 3a, both 
intraband and interband terms are present. The latter 
will be particularly relevant for the susceptibility in the 
presence of spin-orbit coupling (known to be significant 
in the case of Na4lr308 [27]). For the intraband terms, 
q is simply twice kz- Collecting these terms, we arrive at 
the T=0 expression [17] 
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where 



is the effective mass. 
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m the bare mass, and n the band index. Restoring h, Iq is 
equal to ■mgj{ASgs/{4:nhV))'^ [28]. <f>n is a phase angle, 
equal to zero if g is a maximum, 7r/2 if a saddlepoint, 
and TT if a minimum. For the interband terms, we must 
replace Di and D2 by their averages over the two bands. 
For the hyper-kagome case, it turns out that Di and 
D2 are interchanged between the two hole bands. As 



TABLE I: Properties of the hole and electron surfaces of 
Na4lr308 used in calculating the oscillatory coupling, k is 
in units of n/a, m* in units of m, and e in units of ts- 



band 


hp 






Dxy 


m; 


5 


(1,1,1.260) 


-0.971 


-4.182 


-2.465 


0.296 


6 


(1,1,1.260) 


-0.971 


-4.182 


2.465 


0.296 


7 


(0,0,0.312) 


1.522 


2.618 


-0.001 


0.382 



0.6- 



(a) 



0.4- 



hole (R) 



0O.2- 



electron (r) 



0.0- 



0.15 

0.10 

^0.05 

■ 0.00-1 

-0.05 
-0.10 
-0.15 



(b) 



-0.2 0.0 



0.2 0.4 







(A ) 



1 00 200 
Z(A) 



300 



FIG. 3: (Color online) (a) Spinon Fermi surface for Na4lr308 
(a=8.985j4). Spanning vectors are indicated by arrows, (b) 
Calculated oscillatory response from Eq. 12. 



a consequence, the interband contributions in this case 
reduce to the same expression as in Eq. 9, but with 



= 1/lA, 



(11) 



The various parameters are listed in Table I. 

Note that the value of m* for the interband term from 
Table I is 0.239. Summing over the two intraband and 
two interband terms, we find the total hole contribution 
to be twice 0.296 -f 0.239 or 1.07. The ratio of this to the 
electron value of 0.382 is 2.8. Noting that in the present 
case, for all contributions, g is a maximum [29], we then 
find 



I{z) 



-/o^(2.8sin(0.182z)-hsin(0.218z)) 
z^ 



(12) 



where we have absorbed the proportionality constant for 
the mass (which involves the spinon hopping integral tg) 
into Jq, with z (the separation of the two ferromagnets) 
in units of A. This function is plotted in Fig. 3b, where 
we find a periodicity consistent with the hole period of 
34. 5A mentioned above, but with beating clearly present 
due to the electron period. This shows the power of this 
technique to resolve complex Fermi surface topologies. 

We note that although at very low temperatures, 
Na4lr308 exhibits a super-linear specific heat, indicating 
that the spinon surface might develop nodes [12], over 
most of the temperature range, it is consistent with an 
ungapped surface [13]. Therefore, we regard this cubic 
material as a promising one to consider in the context of 
our proposed experiment, though we encourage that all 
candidate spin liquids be looked at. 
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In conclusion, we have proposed an experiment based 
on oscillatory coupling between two ferromagnets with 
a spin liquid spacer for possible detection of a spinon 
Fermi surface. This, and other Friedel like experiments, 
will hopefully be pursued in the future to see whether 
this novel spin liquid state exists. 

Work at Argonne National Laboratory was supported 
by the U.S. DOE, Office of Science, under Contract 
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